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Abstract 

The structural properties of geodesic currents in an ambient Kerr background is studied from 
an analytical point of view. The geodesies in the congruence correspond to charged particles 
that carry energy and angular momentum from the black hole through the Blandford-Znajek 
mechanism. It is shown that the resulting magnetosphere naturally satisfies the Znajek regularity 
condition. Particular attention is paid here to the energy extracted by matter currents rather than 
by electromagnetic Poynting fiuxes. 
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Force-free Electrodynamics of Kerr Black holes 



I. INTRODUCTION 

Nearly thirty years ago, Blandford and Znajek [l| described a method by which a black 
hole could release its rotational energy into the magnetosphere via electromagnetic processes. 
They argued for the existence of a force- free magnetosphere, wherein outflowing currents 
and electromagnetic energy fluxes could in principle carry energy and angular momentum 
away from the black hole. Subsequently, a considerable amount of research has been done 
in trying to understand the nature of the solutions provided by a force-free electrodjTiamics 
in the vicinity of a rotating black hole (e.g., 0, ^, 0, ^). Such solutions are extremely 
important from an astrophysical point of view, as they would be leading candidates for the 
explanation of jets emanating from black holes. Owing to the complexity of the problem, 
recent efforts have been primarily numerical in nature (e.g., 3,0])- 

After obtaining an approximate solution to the Blanford-Znajek mechanism of energy ex- 
traction from rotating black holes that generalizes the original perturbative split monopole 
solution 9], we were recently successful jl^ in providing a class of exact analytic solutions 
to the equations of force-free electrodynamics in Kerr geometry. Although finite everywhere 
in the magnetosphere, these solutions were not, however, physically realistic, but provide a 
basis for further analytic study. For the case of stationary, axisymmetric, force-free electro- 
dynamics in a Kerr background, we 10| showed that that the exact solution permitted an 
electromagnetic current 4- vector of the form 

- - ''jA^] r . (1) 



a?^/^d9 sin 6 

Here, a is the angular momentum per unit mass of the black hole, g is the usual determinant 
of the Kerr spacetime metric g^y, and A is an arbitrary function of 9. Also is the infalling 
principle null geodesies of the Kerr geometry, given explicitly in Boyer-Lindquist coordinates 
{t,r,^,(^} by 

/'^ = (^— ,-1,0,-). (2) 

Here, A = + — 2Mr, where M is the mass of the black hole. Since the geodesic in 
eq. ([2]) is infalling, it is not possible to extract energy from the hole. More importantly, the 
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current flows along null geodesies. While we expect the currents to flow through geodesies 
under force-free conditions, the fact that the geodesies are null is a less desirable feature. 

In this paper, we study the general structural properties of geodesic currents of particles 
with mass when the magnetosphere is force-free. Instead of working with the covariant 
brmalism, we will write all equations of electrodynamics in an intuitive 3 + 1 formalism 
6|. After a brief introduction to the 3 + 1 formalism in Section [Tll we describe the relevant 
equations of a force-free, axisymmetric, stationary magnetosphere of a Kerr black hole in 
Section UTII The structural properties of geodesic currents are described in Section HVl where 
we demonstrate that the Znajek regularity condition is satisfled. Energy and angular mo- 
mentum extraction through particle currents in the force-free magnetosphere are considered 
in Section |Vl and we show in Section |Vl] that a naked singularity cannot result from this 
process. We conclude in Section IVlTl 

II. 3+1 ELECTRODYNAMICS IN KERR GEOMETRY: A PRIMER 

In the 3+1 formalism of electrodynamics, we rewrite the covariant Maxwell equations 
in terms of 3-vectors E, B, D, and H. While these quantities may loose meaning under 
arbitrary coordinate transformations, in a particular frame we have the advantage of working 
with the four familiar equations of electrodynamics. Much of the background and details 
presented in this section can be found in 6] and j^. 

The 3-vectors are deflned on spacelike slices of spacetime. We shall foliate the Kerr 
spacetime by a collection of absolute space described by surfaces of constant time t in the 
Boyer-Lindquist coordinates {t, r, 6, ip} of the Kerr geometry. Here, the metric takes the 
form 

ds'^ = -a^)dt^ + 2 (3^dipdt + 7„rfr^ + 'jeedO^ + -f^^dip^. (3) 
The metric coefficients are given by 

^2 2 . 2Mr ^ -2Mrasm'^e 

(j -a = g^^ = -l^ — , (j = = , = _ , 

A 

2 ^ S^sin^^ 
796 = p , and 7^^ = — — , (4) 

where 

p2 = ^2 ^ ^2 ^Qg2 Q ^ A = - 2Mr + and = (r^ + a^Y - A sin^ 9 . (5) 
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Additionally 

o? = = ^ , and = a ^ = sin ^ . (6) 

Parameters M and a are the mass and angular momentum per unit mass respectively of the 
Kerr black hole. The three-dimensional space for some fixed value of t is endowed with the 
metric 

^^\v ^ IrvAr^ + leeAQ"^ + lif><pdip^. (7) 

Since the spacetime is stationary, the value of t for the absolute space will not matter in 
our discussion. The covariant Maxwell's equations are 

V^F"'^ = , and V^F"'^ = . (8) 

Here is the Maxwell stress tensor, and is the four vector of the electric current. In 

order to rewrite the covariant equations of electrodynamics in a familiar Maxwell-type form, 
three vectors D and H in our absolute space are defined by 

D' = aF^' , and = ^ae^feF^^ (9) 

Here, eijk is the completely anti-symmetric Levi-Civita pseudo-tensor of the absolute space 
such that Sreif — y/j, where, 7 is the determinant of ds^p. In addition, we define the shift 
3-vector as 

P = (0, 0, P^) . (10) 

The constitutive equations relating the electromagnetic field vectors, E and B, and their 
duals, D and H, for regions of zero electric and magnetic susceptibilities are 

E^aD + Px B , and H ^ aB - P X D , (11) 

where the curl is defined by the equation 

{A X By = e'^'^AjBk . 

With these definitions. Maxwell's equations become 

V • S = , and + V X F = (12) 

for the homogenous part and 

V -D^p^, and -dtD + V x H ^ J (13) 
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for the inhomogeneous part. Here, 

Pc = aP , J'' = al^ , (14) 

are the charge and current densities respectively, V is the covariant derivative induced by 
the spatial metric 7jj on our absolute space, and a is the function that appears in eq. (jnD- 



III. STATIONARY, AXISYMMETRIC, FORCE-FREE ELECTRODYNAMICS 

In the Blanford-Znajek mechanism, the magnetosphere is assumed to be force-free. In ad- 
dition, Blandford and Znajek impose, for simplicity, the inherent symmetry of the geometry 
of the rotating black hole onto the currents and fields, namely, stationarity and axisymmetry. 
Consequently the Lorentz force on the currents in the magnetosphere is trivial, so that 

F^J" = . (15) 

In our 3 + 1 language, the above condition takes the form 

E-J = 0, (16) 

and 

p^E + JxB = 0. (17) 

The consequences of eqs. ( ITBi) and (fT7|) along with the requirements of stationarity and 
axisymmetry will be summarized below. The interested reader is referred to [gl] for further 
details. The poloidal and toroidal components {Ap and At, respectively) of a vector field 
are defined such that A = Ap + At, where Ap = A^dr + A^de and At = A^d^. Since the 
magnetosphere is stationary and axisymmetric, 

= . (18) 

This is true because stationarity implies that E is the gradient of a scalar, and for axisym- 
metric fields, all (p derivatives vanish. From eq. (IT7|) we have that 

E ■ 5 = 0. (19) 

Eqs. f fTSj) and f|T9|) imply that there exists a vector oj = VL d^p such that 

E = -ujxB. (20) 
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Here, is a function dependent on r and 6. Since B is divergence free, 

B = V X A, (21) 

where A is the 3-vector potential associated with B. Therefore, 

Bp = ^{A^^edr - A^^rde) , (22) 

since A^p is independent of t and (p. Surfaces of constant A^ are called equipotential surfaces. 
Here is the toroidal component of the 3-vector potential A. Functions that are constant 
on equipotential surfaces are called equipotential surface functions, i.e., a function / is an 
equipotential surface function if 

Bpf = = A^^g drf - A^,. def . (23) 

It can be shown that VL and H^p are equipotential surface functions, and consequently from 
eq. ( l20i) we get 

Ep = -d [n dA^ . (24) 



Since Q is an equipotential surface function, it is strictly a function of A^, and hence the 
above integral is well defined. The constitutive eqs. (|TT1) give 

D = Dp = -^{Q + P'^)dA^, (25) 



and 



Hp={a'-I3^-I3^Q)^ . (26) 



The electric charge density is determined by the divergence of Dp: 

- \flPc = <9r[^^(7w^ + Pv) lee A^^r] + <94^^(7w^ + l^v) 7rr A^fi]. (27) 



The toroidal current vector can be obtained by computing the component of the curl in 
eq. ([T3D: 



(28) 

The poloidal current is completely determined by H^p and is given by 

y/yjp = H^fidr - H^^rde- (29) 
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Finally, 

= H^/a. 

Since is an equipotential surface function, from eq. ( l29l) we have that 



Jp = -^^{A^,edr - A^^rde) = , (30) 

dA^ dA^ 

i.e., the currents flow along equipotential surfaces. The only remaining requirement in the 

explicit formulae given above for the fields and current stems from eq. (fTTl) . and can be 

written in the form 

= -«(Pc^^7w-^) • (31) 
Since the left hand side of the equation above is an equipotential surface function, we must 
pick an equipotential surface function Q such that the right hand side of the equation is an 
equipotential surface function as well. Additionally, for the fields and currents to remain 
well defined at the event horizon of the Kerr black hole, given by 



r = r+ = M + - a2 



we must have that 

\^^= !^5'^(2Mr n-a) |,+ (32) 



a 



in Boyer-Lindquist coordinates. This is the Znajek regularity condition derived in Hj]. In 
the above equation, all quantities are to be evaluated at the event horizon. 

The poloidal component of the Poynting vector for the electromagnetic field is given by 

Sp = {Ex H)p . 

Consequently, the rate of extraction of energy from the black hole is given by 

dE 

p.m. 



dt 



j nn^B'^ de dcp . (33) 

J r9= 00 



' r2= 00 

In a similar manner, using the fact that is a Killing vector field of the geometry, we get 
the following expression for the rate of extraction of angular momentum from the black hole: 

dLp 



dt 



I H^B' ^dOdip. (34) 



' r2= 00 

Under this circumstance, the mass and angular momentum of the black hole decreases at a 
rate given by 

5M dEem , 5J dLem 

—r- = ; — , and — = — . (35) 

5t dt 5t dt ^ ^ 
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IV. THE GENERAL FORM OF GEODESIC CURRENTS 



In the case of force-free electrodynamics in a curved spacetime, we would expect the 
charged particles to flow along the geodesies of the geometry (a fact that was manifest in 
eq. ([I])). With this in mind, we write the current vector in the form 

r = F{r, 6) u" . (36) 

Here, u'^ is a timelike geodesic of the Kerr geometry. The function F relates the 4-velocity 
of the charged particle at a point to electromagnetic current density at that point. The 
properties of u" are well understood (for example see [l^ and 13]). In particular, u" can 
be written in the form 

= {t,^,^,^) . 37 

Here, t, 0, a/R and are functions of r and 6. Also, in ^/R and we haven't yet com- 
mitted to the "+" root. The notion of time (dot) derivative is not relevant to our analysis; 
the notation has been simply borrowed from the existing literature on Kerr geodesies. The 
explicit forms of the geodesic functions are 

. E'^S-2aMrC 



and 



_ 2aMr£ sin^ + {p^ - 2Mr)C 

^ " JK^e ■ ^^^^ 



The functions R and are given by 

R = C^ + A{q\^ - /C) (40) 

and 

e = /C + q^a^ cos^ 9 ^ , (41) 

sin 6 

where 

C = (r^ + a^)S -aC and D = C - aS sin^ 6 . (42) 
In the study of Kerr geodesies, S, C, and /C are simply constants for a given geodesic, and 

= 1, 0, -1 
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for spacelike, null, and timelike geodesies respectively. The terms 8 and C can be related to 
the energy and angular momentum of the particle that flows along the geodesic. Explicitly, 



8 = -g{K, u) = -{gut + gt^'p) (43) 

and 

C = g{m, u) = {gt^i + 7<^^0) . (44) 

Here K = dt and m = d^p. K, is the famous Carter's constant [l^. In our case, since 
the currents flow along equipotential surfaces, 8, C and /C become equipotential surface 
functions, which are now subject to the condition 

/,r + /,e = , (45) 

for any equipotential surface function /. Also, since we are only interested in particles 
with mass, we will set = —1. Conservation of electric charge is implied by the equation 
V^/^ = 0, i.e.. 



dr[F^y^+de[F^V^ = . (46) 

This is actually an equation for F. 

In writing eq. fl36p . we have assumed the explicit form of the charge density and the 
3-vector current density. In the remainder of this section, we will construct a formula for 
H^p and show that it satisfies the Znajek regularity condition, eq. fl32|) . But first we present 
a few preliminaries. 

From eq. (|T^ we see that 



Pc = OiF i , J"^ = aF ip , 

and 



Jp = aF + aFy^de. (47) 



Consequently, eq. (]3Ti) can be written as 



1 dH"^ 

-^ = Asm'eF{in-p). (48) 



In the theorem that follows, we obtain an explicit expression for H, 
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Theorem 1. In force-free, stationary, axisymmetric electrodynamics in a Kerr background 
when the electromagnetic current vector takes the form given by eq. fl36|) . 

^^^A=m«JU)-^^^_^ (49) 

Proof: From eqs. (^U^ and fH7|) we get that 

H^ ^ = -F sin eVe . (50) 

But note also that 

TT TT f_ A 

Substituting eqs. (148|) and (!50|) into the above equation, we get the necessary result. □ 

Since is an equipotential surface function, if^ can also be written in the form 

H, = zA^R^p^A^^,, (51) 
V R 



Corollary 1. as given by eq. (149|1 and/or eq. (I5T1) will be an equipotential surface 
function if and only if 

Rdri^) + edei^)=0, (52) 
V R V H 



where 

X = AsinO {i n - . 

Proof: From eqs. (H9!) and (15T|) . 

rr 4 4 

and since if^ is to be an equipotential surface function 

y/R H^^r + VeH^,e = 0. 

Since is an equipotential function, the above two equations immediately imply eq. (l52l) . 

□ 

One of the main advantages, and justification for eq. fl5Bl) is that the expression for if^ 
as derived above naturally satisfies the Znajek regularity condition. 
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Theorem 2. The form of H^p as given in eq. (15T!) (or equivalently in eq. (H9!) ) satisfies the 
Znajek regularity condition, eq. fl32|) . when 

i?L , = — C . 



Note, the only choices for a/R|,.^ are ± C. 
Proof: From eq. ([22]), 



Also, 



_ - /\sme {t Vt - Lp) 



We choose the square root such that 

V^|a=o = — C|a=o • 
Here, C is the function defined in eq. fH2|) . Also, 

p^A i|A=o = (r+ + a^f£ - 2aMr+ C = 2Mr+ C|a=o , 

and 

p^A 0|a=o = a C|a=o • 

Putting all the above equations together, we get the stated result. □ 

V. THE EXTRACTION OF ENERGY AND ANGULAR MOMENTUM VIA PAR- 
TICLE JETS 

In Section IllH the terms representing the extraction of energy and angular momentum 
(eqs. (l33il and were derived from the conserved quantities of the electromagnetic stress 
tensor. In this section, we will do the same for the matter fields. This will be useful 
for calculating the energy released via matter jets from supermassive black holes. Separate 
conservation of field and particle energy and angular momentum stems from the fact that we 
are only considering a force-free magnetosphere. In particular, the total energy-momentum 
tensor of magnetohydrodynamics can be written in the form: 



rpfiu rp^u I rp/liy 



matter 
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Since this must be divergence free, we must have that 



But, from eq. f|T5|) we see that 



Consequently, we must have that 



V T^^" = 



V T^'^ — n 

^ fj.-^ matter ^ ) 



as mentioned above. When the magnetosphere is filled with a geodesic congruence, we can 
take the matter energy momentum tensor as 

TZtter = (m/e) Fir, 9) , 

where e and m are the electric charge and mass of the particle species, respectively, and F 
and are as given in eq. (136|) . For then, 

* M matter 

(m/e) V/,(F u^)u'' + (m/e) F u^^VX = • 

In the equation above, the first term on the right hand side vanishes due to the continuity 
equation for the electromagnetic current density, and the second term is trivial because of 
the geodesic equation. If the charge of the currents are position dependent, we must either 
take the derivative of the e/m term in the above equation, or consider the regions of different 
values for e/m separately. 

A. Matter Energy Flux 

Since K is a. Killing vector field of the Kerr geometry, 

i'^u matter ^ ) ~ ^ i 



V / I mailers • r— 



Let 



e = -a T/^„,jg^ (53) 
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be the energy density of the matter fields, and 

S^^-a Ti^^^ (54) 



be the matter energy density flux. Therefore, 

[edV^-[v-SdV^-[ S-ndA. (55) 

Jv Jv JdV 



d 

di 



Here, y is a 3-dimensional manifold in our absolute space, dV is the boundary of V, and n 
is the unit outward pointing normal on dV as per the divergence theorem. Consequently, 
the rate of matter extraction of energy from the black hole is given by 



dEm 

dt .1 



I y/^dOdcp. (56) 

J ri= oo 

Clearly, 

Sr^-a-F^{gui + gt^^). (57) 
e 

But 

-{gtti + gt^ (p) 
is nothing more than the poloidal function £, therefore. 



dEjn 



= — / F VR£ sin 9 de dip . (58) 

^ J ri= oo 



dt 

B. Matter Angular Momentum Flux 

Using the Killing vector fleld m = d^, we can deduce a similar expression for the rate of 
extraction of angular momentum from the black hole via particle jets. Setting 

l^(^Tl (59) 

as the angular momentum density of the matter fields, and 

L'^^^T^ matter (60) 

as the matter angular momentum density flux, we see that the rate of extraction of angular 
momentum from the black hole via particle jets is given by 



dt e 



j F \fRC smO de d^ . (61) 

J ri= oo 
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VI. THE AREA THEOREM, AND THE BLANDFORD-ZNAJEK MECHANISM 

In it was shown that for electromagnetic extraction of energy from the horizon of the 
black hole, the rate of energy extraction was related to the rate of angular momentum 
extraction by the expression 

Here, 

is the angular velocity of the event horizon. A black hole that is adiabatically evolving along 
a Kerr sequence such that its rate of energy and angular momentum extraction satisfies 

dE ^ dL 



cannot lead to the formation of a naked singularity 15(]. This can be easily seen as follows. 



The irreducible mass Mj^r of a Kerr Black hole is given by 



Ml^ = -[M' + VM^-J^], (65) 

where J = aM is the angular momentum of the black hole. Under extraction rates subject 
to eq. flM|) . the mass and the angular momentum of the Black hole decreases to 

M ^ M + 5M and J ^ J + 6J , (66) 

where SM and SJ satisfy 

6M >nH6J . (67) 
We now state the following calculation due to Christodoulou [15]) without proof. 

Theorem 3. When energy and angular momentum is extracted from a Kerr black hole so 
that eq. (|671) (or equivalently eq. (16^ ) is satisfied, we have that SM^.^ > . 

The area of the event horizon is given by 

^ = / y/geeg^ dOdip = IGvrMf^^ . (68) 

J r-i- 

Since the area of the event horizon is proportional to Mf^j., from the theorem above, we see 
that it is non-decreasing, and so, we have the following corollary. 
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Corollary 2. A black hole that is adiabatically evolving along a Kerr sequence such that 
its rate of energy and angular momentum extraction satisfies eq. fl64l) always contains an 
event horizon that cloaks the physical singularity of the Kerr geometry. 

Of course, when all the angular momentum is extracted, the Kerr black hole reduces to 
the Schwarzschild solution, and here too we do not have a naked singularity. A dynamically 
evolving black hole, in general, need not follow a Kerr sequence. In this case, the powerful 
area theorem of Stephan Hawking requires that the area of the event horizon continues 
to be non- decreasing. 

We conclude our analysis of the geodesies currents in the force-free magnetosphere of the 
Kerr geometry by showing that the matter energy and angular momentum extraction rates 
are subject to the same inequality as their electromagnetic counterparts (eq. flMl) ). and 
hence we are in no danger of violating the Penrose conjecture. 

Theorem 4. When the currents in a force- free magnetosphere are given by eq. fl36l) . the 
extraction of energy from the matter fields will not lead to the formation of a naked singu- 
larity. 

Proof: At the horizon, if charge are to fiow out to the magnetosphere, from theorem ([2]) 
we need that 

y/R\^^ = -C>0 . 

This is possible only when the geodesies are such that 

8 <VLh (69) 
and consequently from eqs. dSSD and fl^ 

Therefore, corollary ([2]) implies the necessary result. □ 



VII. CONCLUSIONS 

The mechanism to form jets from black holes is a major unsolved problems in astrophysics, 
and its explanation is important for our understanding of nonthermal radiation form Solar- 
mass black holes in the Milky Way, newly born black holes associated with gamma-ray 
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bursts, and supermassive black holes at the centers of distant galaxies. One possibility 
is that black holes with radio-emitting jets are distinguished from radio-quiet black holes 
by extracting rotational energy from the black hole through the Blandford-Znajek process. 
Matter currents can also carry energy and angular momentum, and represent a plausible 
alternative to electromagnetic Poynting flux for extracting the black hole's rotational energy. 

In this paper, we have used the 3 + 1 formulation of general relativity to analyze currents 
that flow in the force-free magnetosphere near a Kerr black hole. We have rewritten the 
constraint equation for a force-free magnetosphere, eq. (13T|) . to apply to particles that follow 
timelike geodesies, yielding eq. ( HHl) . The toroidal magnetic susceptibility, if^, from which 
the poloidal current is derived, is shown to satisfy the Znajek regularity condition. Ex- 
pressions for energy and angular momentum flux associated with the matter currents were 
obtained. We show that extraction of energy and angular momentum cannot lead to the 
formation of a naked singularity. In future work, trajectories of particles in the black hole 
magnetosphere will be derived, and will be compared with observations of jets from black 
holes. 
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